Particle-particle random phase approximation applied to Beryllium isotopes 
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This work is dedicated to the study of even-even ~^ Be isotopes using the particle-particle Ran- 
dom Phase Approximation that accounts for two-body correlations in the core nucleus. A better 
description of energies and two-particle amplitudes is obtained in comparison with models assuming 
a neutron closed-shell (or subshell) core. A Wood-Saxon potential corrected by a phenomenological 
particle-vibration coupling term has been used for the neutron-core interaction and the DIS Gogny 
force for the neutron-neutron interaction. Calculated ground state properties as well as excited 
state ones are discussed and compared to experimental data. In particular, results suggest the same 
2si/2-lpi/2 shell inversion in ^'^Be as in ^^Be. 



I. INTRODUCTION 

In two previous studies [l|, y] j it has been proposed the 
same inversion in ^'^Be between 2si/2 and lpi/2 shells 
as the one found in ^^Be Q and ^°Li Q. This assump- 
tion has been suggested to obtain a good description of 
^^Be two-neutron separation energy. In Ref. l| , a simple 
pairing model was utilized whereas in Ref. [2] a particle- 
particle RPA (pp-RPA) approach with the DIS Gogny 
force [5|, Q was introduced. 

The present study is an extension of the work of 
Pacheco and Vinh Mau [2]. It aims at getting addi- 
tional and essential information on the structure of ^'^Be 
from the study of ^^Be and ^^Be. Recently, several mea- 
surements on Beryllium isotopes and especially on their 
excited states have been performed [7H13|. These new 
experimental data associated with our present calcula- 
tions are important both to assess the validity of the 
model itself for the description of excited states and to 
test with more constraints the hypothesis of the inversion 
between 2si/2 and lpi/2 shells in ^'^Be. The pp-RPA ap- 
proach is a three-body model that provides information 
on {A — 2), A and {A + 2) nuclei starting from a reference 
core nucleus with A nucleons. A spherical symmetry is 
assumed. In the present work, we apply it both to ^"^Bc 
and ^^Be cores in order to get consistent information on 
^"^"^Be isotopes. We use a Wood-Saxon (WS) poten- 
tial corrected by a phenomenological particle-vibration 
coupling [l|, y, [iJl for the neutron-core interaction and 
the DIS Gogny force to describe the neutron-neutron ef- 
fective interaction. Moreover, in the previous study of 
Pacheco and Vinh Mau Q], the two-body spin-orbit in- 
teraction term was neglected in the study of 0^ ground 
state [la]- As in the present work we are interested in 
excited states with different spins and parities, the spin- 
orbit term is included. 

In this article, results associated with ground state 
properties of even-even *~^^Be isotopes as well as their 
excited spectra and transition probabilities are presented. 
Comparisons with experimental data are also discussed 



in detail. In Sec|TTl we briefly introduce the pp-RPA 
approach and give a few analytical formulas concerning 
radii and transition probabilities. Sec lIIII is devoted to 
results obtained from a ^"Be core. The neutron-^°Be in- 
teraction has been fitted from experiment [l6| and used 
to generate the single particle basis. The ability of the 
model to reproduce the experimental knowledge on ^"^Be 
and ^^Be is demonstrated. In Scc lIV| similar analysis 
with a ^■^Be core is presented. As the neutron-^^Be in- 
teraction is not known precisely from experiment, we 
investigate different scenarii for the neutron-core inter- 
action, constrained by experimental knowledge on ^^Be, 
^■^Be and ^'^Bc. Conclusions are given in SecfVl 



II. PP-RPA FORMALISM AND ASSOCIATED 
OBSERVABLES 



The pp-RPA is a well-known formalism used to study 
nuclei which can be approximated as a core plus or mi- 
nus two nucleons. An important property of this ap- 
proach is its ability to account for two-body correlations 
in the core nucleus. Different ways of deriving pp-RPA 
equations can be found in the literature [l7|. For ex- 
ample, the Green's function method shows that pp-RPA 
approach introduces in the core ground state contribu- 
tion of multiparticle-multihole configurations [l8| . 

In this part, we recall briefly the standard equations in 
order to introduce our notations. Starting from a nucleus 
with A nucleons, the pp-RPA equations describing the 
nuclei with {A + 2) and {A — 2) nucleons are. 
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{n - [e„ + e/j]) Xap + y^^{kl\V\al3)xki 

kl 

+ "^{k\\V\^P)x^x = 0. (2) 

In Eqs.(IT]) and ([2]), the set of e arc single particle en- 
ergies determined together with the corresponding wave 
functions in a given neutron-core potential. The cho- 
sen potential is presented in Sec lIIII (see Eq. pT]) ). Ma- 
trix elements of the two-body interaction V are anti- 
symmetrized. Latin indices stand for unoccupied single 
particle orbits and greek indices for occupied ones. The 
quantities x gives the standard two-nucleon amplitudes 
X and Y for nuclei with {A + 2) and {A — 2) nucleons, 
respectively. They read 
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where m <n. The states of {A + 2) or [A — 2) nuclei are 
labeled by N or M, respectively. The state |A,0) stands 
for the correlated ground state of the core. \A + 2,N) 
and \A - 2, M) are the states of the {A + 2) and (A - 2) 
nuclei, respectively. In Eq.Q, the pair creation operator 
A'^ is defined as 

■Kibi'J'^'^j) ="^ab ^ {jajbmamb\J,Mj)alal, (5) 
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with the two nucle- 



with a < b and Vab = (1 + ^jajb) 
ons in the same spherical j-orbital. The annihilation op- 
erator A in Eq. @ is deduced from Eq. (IS|) . These expres- 
sions are valid for both occupied and unoccupied states. 
In Eqs.dH and ([2]), il is the energy of the state related 
to {A + 2) or {A — 2) nucleus taking as reference the 
A nucleus ground state energy. Hence, for the {A + 2) 
nucleus. 
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and for the {A — 2) nucleus. 
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More details concerning the pp-RPA formalism are given 
in Appendix VK\ 

In addition to energies and amplitudes, other observables 
as rms radii and transition probabilities will be discussed 
in Sec mil and II VI We give here the expressions for those 
two quantities. Details can be found in Appendices iBl 
andlCl Following Ref . [l3| . the rms radius (r^)^^_2 of the 
{A + 2) system can be expressed in terms of the radius 
of the core (r'^) 
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with 
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In Ea. ([TT|) . A is the distance between the center-of-mass 
of the two extra-nucleons and the center-of-mass of the 
core and p the distance between the two nucleons, 



A = -in 

2^ 
P = ri ra. 
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where ri and r2 are the coordinates of the two extra- 
nucleons relative to the center-of-mass of the core. 

The model assumes an inert core plus two correlated 
neutrons. Therefore the B{E1) for a transition from the 
of ground state to the 17 excited state is given by the 
so-called soft dipole strength [20], 
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(14) 

where the sum runs over the two extra-neutrons. The 
contribution to the El strength comes from only the two 
extra neutrons with an effective charge, 

(15) 



" A + 2' 

where Z is the number of protons in the core. The expres- 
sion of B{E1) and more general transition probabilities 
in the pp-RPA formalism are given in Appendix [Cj In 
addition, a simple expression for the sum of B{E1) over 
all the dipole states can be derived |22l |. 
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This formula is very useful because it provides a con- 
straint between ^^ Bf{El) and A: the El strength ex- 
tracted from experiment should not exceed the value ob- 
tained in the right hand side of Eq. ([T5)) , calculated with 
the experimental value of (A"). 



III. 



DESCRIPTION OF EVEN-EVEN ^-^^BE 
FROM A ^"BE CORE 



In this part, we study even-even ^~^^Be isotopes us- 
ing the pp-RPA formalism with a ^°Be core. It is well- 
established that the pure Hartree-Fock (HP) approxima- 
tion fails to reproduce ^^Be properties, in particular the 
inversion between 1/2+ and l/2~ states [lal- This inver- 
sion is due to the coupling with a low energy 2"*" collec- 
tive state of the -^°Be core (E2;=3.36 MeV) characterized 
by a strong quadrupole transition probability {B{E2;0f 
^■2+)=52 e^fm") O, [HI, HI] . A phenomcnological cor- 
rection to the one-body potential, simulating the cou- 
pling between the neutron and a phonon of the core, has 



been proposed in Ref.|l^ 
is written as 



This new one-body potential 
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In Eq. pTj) . Vws is a WS potential including central and 
spin-orbit parts plus a symmetry energy term accounting 
for neutron excess. In Eg.p^. f{r) is the Fermi form fac- 
tor of the WS potential characterized by the diffuseness 
a. The values of parameters of Ref. [2| have been adopted. 
Eigen wave functions of this adopted WS potential con- 
tain the effect of the coupling between a neutron and a 
phonon of the core. They are no longer pure one-particle 
wave functions but a mixing of one single neutron state 
coupled to the ground state of the core and of a single 
neutron state coupled to a phonon of the core. In prac- 
tice, parameters a^ of Eq. p^ have been fitted from the 
experimental spectrum of ^^Be [l6[ for 2sx/2, lpi/2 and 
ld5/2 neutron states and taken equal to zero for the other 
states. From a technical viewpoint, a 20 fm range box 
has been used to determine eigen wave functions. 

In Table m results obtained for ground state properties 
of ^°'^^Bc are presented. A good agreement for two neu- 
tron separation energies S2n in ^°Be and ^^Be is found. 
The radius of ^^Be calculated using Eg. pH)) is also well 
reproduced. Predictions for the mean distances (A^)^/^ 
and (p^)^/^ associated with A and p, Egs. p^ and ([T^. 
are indicated although no experimental values are yet 
available. 



TABLE II: Main pp-RPA amplitudes of 0^ ground state in 
i^Be. 



Xab 
(lPl/2)' (lPl/2, 2pi/2) (2Si/2)' (ld5/2)' 



0.76 



0.31 



0.50 



0.43 



X, 
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(IP: 



3/2; 



0.57 



experimental data [26[ and is consistent with results pro- 
vided by other models |21| . 

As discussed, results obtained for ground state prop- 
erties of ^"Be and ^^Be are in an overall good agreement 
with known experimental data. Now we extend our cal- 
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FIG. 1: Low lying spectra of Be obtained with pp-RPA 
compared to experiment. 



TABLE I: Theoretical and experimental values of 5*2,1 (MeV) 
in ^°Be and ^^Be (from Ref.lH), (r^^^ (fm), {X^/^ (fm) 
and (p^>y+2 (fm) in ^^Be (from Ref.[2l). 



S2„("Be) S2n(^'Be) 


{rrii, ip^^' i^r^' 


Theory 8.49 3.62 


2.76 4.89 4.10 


Exp. 8.48 3.67±0.01 


2.59±0.06 



In Table |TT1 the main pp-RPA amplitudes contribut- 
ing to the ground state wave function of ^^Be are shown. 
From the values of Xab that describe the contribution of 
a configuration where two neutrons are created in two 
unoccupied particle states, one sees that the wave func- 
tion of ^^Be is a mixing of different configurations with 
the two extra-neutrons mainly in (lpi/2)^, (2si/2)^ and 
(ld5/2)^ configurations. The value of Aq,^ that stands for 
two-neutron configurations in hole states is quite large 
indicating strong correlations in the ^"Be-core with the 
presence of two particles - two holes configurations. Such 
a strong mixing is an indication of a breakdown of the 
N=8 shell closure in -'^^Be. This result is supported by 



culation to the description of excited states. We focus on 
the 0+, 0~, 1~ and 2+ excited states of *Be and ^^Be. 
Calculated excitation energies E^; as well as dipole tran- 
sition probabilities B(E1) are compared to experimental 
ones. In practice, the pp-RPA formalism provides the 
energy E^ of a state relative to the two neutrons + core 
threshold and the two neutron separation energy S2n. 
The excitation energy E^; is thus expressed as 



Err — Er 



S2 



(19) 



Concerning ^Be, our model predicts the absence of low 
lying 1~ and 0+ excited states that is in agreement with 
experimental data. Only a 2+ state with an excitation 
energy of 3.82 MeV is found. This state is experimentally 
observed at a slightly lower energy of 3.03 MeV. 

Results for ^^Be are summarized in FigJT] The experi- 
mental Oj excited state at E2:=2.24 MeV [3| is well repro- 
duced with a theoretical excitation energy of 2.48 MeV. 
In Table IIII[ the main amplitudes corresponding to this 
state which is mainly built from (2si/2)^ and (lpi/2)^ 
configurations are shown. It can be noted that for O2 



TABLE III: Main pp-RPA amplitudes of Oj state in ^^Be. 
X(i3 amplitudes are found negligible. 



-0.48 



Xa(, 
(2Si/2)^ 

-0.86 



(lPl/2,2pi/2) (2Si/2,3Si/2) 

0.11 0.10 



also recall that in -'^^Li the calculated 17 low lyins; state 
[29| was in agreement with later measurement [30|, both 
for excitation energy and transition probability. In ad- 
dition, our model predicts two higher 1^ states charac- 
terized by E^=4.24 MeV with S(£'l)=0.064 e^ fm^ and 
E:r==4.32 MeV with B(£'l)=0.066 e^ fm^. These states 
have not been yet observed experimentally. 



state no effect of correlations in the core, which are repre- 
sented by negligible contributions from Xq/j amplitudes, 
are found. Looking at both Tables |ll] and IIIH one sees 
that the model gives a higher contribution of the (2si/2)^ 
configuration in the O2 excited state than in the ground 
state 0^ . This trend may be related to the results of Ka- 
nungo et al. [13| where the s-wave spectroscopic factor 
was equal to 0.28 for the ground state and 0.73 for the 
0^ excited state. 

We also obtain a lj~ state at Ea;=2.59 MeV with a tran- 
sition probability of B{El;Ot -^ ir)=0.45 e^ fm^. The 
associated experimental values are Ea;=2.68(3) MeV and 
B{El)=0mi{13) c2 fm^ [23. The energy of Ij" state 
is in agreement with experiment whereas the value of 
the calculated B{E\) is overestimated by a factor of 10. 
Sagawa et al. [28| find a low energy transition strength of 
i3(£'l)=0.063 e^ fm^ in the context of large shell model 
calculations using extended wave functions for loosely- 
bound lpi/2 ans 2si/2 states. In our model, as equations 
are solved in a box, we are free from this kind of cor- 
rection as our cigen wave functions have already good 
asymptotic properties. One difference may come from 
the fact that Sagawa et al. renormalize the depth of 
their HF potential to reproduce half of the empirical two- 
neutron separation energy. Their deduced single particle 
wave functions are thus much more bound than ours, 
with certainly a smaller spatial extension. This may ex- 
plain partly why they obtain a lower El strength. 

In Figl^l the calculated El strength distribution in 
^^Be is shown. The El strength associated with the 1^ 
state gives the largest contribution in our calculation. In 
Ref.[2^ where no core is assumed, the contribution of 
the giant dipolc resonance is found in the energy range 
E2:=10-13 MeV. In our calculation only the soft dipolc 
part of the El strength is accessible. With the help of 
Ea.([T51) and the calculated value of A, the total deduced 
El strength is equal to 1.8 e^ fm^. Transition proba- 
bilities B{E1) as well as radii are two types of observ- 
ables very sensitive to the content of the wave function. 
For ^^Be, the ground state is a large mixture of (2si/2)^, 
(lpi/2)^ and (ld5/2)^ configurations, as shown in Table 
im while the 1~ states are nearly pure two-neutron con- 
figurations. In particular, the 1^ state is nearly a pure 
(lpi/2, 2S1/2) configuration. As discussed earlier, the 
calculated ground state wave function of ^^Be provides a 
good value for the radius which depends strongly on the 
values of A and p. This agreement strongly lends cre- 
dence to the calculated value of J2f ^fi^^)- One may 
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FIG. 2: Calculated results of El strength distribution in ^^Be. 



Regarding 2+ states, our model is unable to repro- 
duce the low lying 2'^ state, experimentally found at 
2.11 MeV i31|. However, it predicts two additional 2+ 
states at higher excitation energies, 3.86 MeV and 4.59 
MeV. These results are coherent with those obtained by 
Romero- Redondo et al. [32|, [S^ where they studied ^^Be 
within a three-body model with an inert core. They were 
unable to reproduce the low 2^ state and also predicted a 
higher state close to our state at 3.8 MeV. Indeed in our 
model the phenomenological particle-vibration coupling 
in Eq. p^ takes into account the effect of the 2^ state of 
^"Bc on the single particle scheme of ^^Be. But the 2f 
state of ^'^Be itself is not present. As shown by Nunes et 
al. [2l|, in order to explain this state, one has to take 
into account explicitly the excitations of the ^°Be core. 

Concerning the Oj" state, it is obtained in our study 
with a relative energy of -0.71 MeV (E3;=2.91 MeV) close 
to the results of Romero-Redondo et al. [321 (Ea;=2.5 
MeV). It is built from a nearly pure (lpi/2, 2S1/2) con- 
figuration. This Oi state has not been yet observed even 
in a recent experiment jl3| . 

As discussed previously, results obtained for ground 
states and excited states are found to be in a quite good 
agreement with experiment, except for the dipole transi- 
tion probabilities and the 2^ state energy. The pp-RPA 
approach is able to describe quite well the known spec- 
trum of *Be and ^^Be and suggests the presence of higher 
states which have not been yet observed experimentally. 



IV. DESCRIPTION OF EVEN-EVEN """BE 
FROM A ^^BE CORE 



In this part, we are interested in the description of 
Be isotopes using the pp-RPA approach 
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starting from a ^^Be core. Concerning neutron states in 
the field of ^^Be, the situation is still not clear. A 5/2+ 
resonance has been first observed at 2.0 MeV above the 
neutron-^^Be threshold [3J]. This state is interpreted as 
one neutron in the ld5/2 shell. Later experiments have 
confirincd this resonance with a relative energy between 
2.0 McV and 2.4 MeV [i-[lS[il,[Mll3- A lower state has 
been observed close to 0.3 MeV [33 • A 1/2+ assignment, 
as suggested in Ref.[37|, with a pure neutron s-state im- 
plies that the shell order in ^^Be is given by a simple WS 
potential with an occupied lpi/2 shell. In a recent exper- 
iment by Kondo et al. [IJ], they have identified -'^■^Be res- 
onances of 1/2" with £'^=0.45 MeV, 1/2+ with £,.=1.17 
MeV, and 5/2+ with £,.=2.34 MeV. This experimental 
result suggests an inversion between the 2si/2 and lpi/2 
shells as the one predicted in Ref.Q]. In order to clar- 
ify those contradictory results, in the following we test 
two sccnarii in ^■^Bc, concerning shell ordering in order 
to find a scenario that reproduces at best experimental 
observables for ^*Be. 

In the first scenario, we assume a normal order of shells 
with a low lying 2si/2 neutron state. The lpi/2 shell 
is then the last occupied neutron orbital in ^^Be with 
an energy e(lpx/2)=-3.17 MeV given by the measured 
neutron separation energy in ^^Be [38|. In the following 
this scenario is referred as scenario A. 

In the second scenario, the inversion between 2si/2 and 
lpi/2 shells is assumed. Indeed, in ^^Be a 2^ state with 
an excitation energy of Ea:(2+)=2.6 MeV and a transi- 
tion probability B(E2;0f -^ 2+)«50 e^ fm" [H is ob- 
served close to the one in ^°Be [ijj 123]. A similar ef- 
fect on the neutron states in ^^Be as the existing one 
in ^^Be can be expected. Thus assuming that the shell 
inversion present in ^^Be holds in ^"^Be, the last occu- 
pied nucleon orbital in ^^Bc is the 2si/2 shell with an 
energy e(2si/2)==-3.17 McV. The lpi/2 shell is unbound 
with an energy not established experimentally. In this 
second scenario the energy of the lpi/2 shell is consid- 
ered as a parameter. The energy of the ld5/2 state is 
assumed to be e(ld5/2)=2.27 MeV, a bit more than the 
usual one (e(ld5/2)=2.0 MeV) but still in agreement with 
recent experiments [10, Il2l ■ This result is not in contra- 
diction with experimental knowledge. Indeed, it has been 
shown in Rcfs.[l3, |40l that a lpi/2 state above threshold 
is needed in order to reproduce the experimental neutron- 
^^Be spectrum. In the following, this second scenario is 
referred as scenario B. 

Now, we compare results obtained with the two sce- 
narii with experimental data on even-even ^''"^^Be, both 
for ground state and excited states. 

In scenario A, we first assign a relative energy of 0.3 
MeV to the 2si/2 state, according to the experimental 



TABLE IV: Theoretical and experimental values of S2n 
(MeV) in ^^Be and "Be (from Ref.lH), {v^)]ll^ (fm), {X^)^''^ 
(fm) and (p^)^''^ (fm) in "Be (from Refs. ^^, A deduced 
using Eq.Uni) in the two cases of non-inversion (A) and in- 
version (B) of the 2si/2 and lpi/2 shells. 





S2„(^'Be) 


52„("Be) (r^>Y^\ (p^)^/^ (A^)V2 


A 
B 


2.91 
3.71 


0.51 3.45 8.45 5.45 
1.29 2.91 4.56 4.02 


Exp. 


3.673±0.015 


1.26 ± 0.01 3.10±0.15 5.4±1.0 4.2±1.7 



suggestion of a low 1/2+ neutron state from Ref.[33|- In 
this case all calculated quantities disagree with experi- 
mental values. In particular, ^^Be is found under-bound. 
We then decrease the energy of the 2si/2 state preserv- 
ing the agreement with experimental data. The energy of 
the ld5/2 state is fixed at 2 McV and the one of the 2si/2 
state at 0.09 MeV (as low as possible ensuring also an un- 
bound ^■^Be). Results are summarized in Table ITVl One 
sees that, even in that case, it is impossible to describe 
correctly both the two-neutron separation energies S2n of 
^^Be and ^"^Be. The rms radius is overestimated and the 
rms value of A is in the upper part of the experimental 
error bars. The rms value of p is also overestimated by 
more than 2 fm. It is interesting to note that if the 2si/2 
state is bound, results are improved and agree with the 
'S'2n(^^Bc) value given by Descouvcmont et al. (43 - l44| . 
If the 2si/2 state is unbound and the energy of the d5/2 
shell is decreased, results closer to experiment arc found, 
as in the work of Thompson and Zhukov [45| . However, 
these assumptions are not justified since a bound 2si/2 
state and a ld5/2 state below 2 MeV disagree with all 
experimental measurements. In addition, scenario A has 
also been studied within a model introducing a core de- 
formation [46| . Only for a very high deformation param- 
eter (/3 > 0.8), ^^Be has an unbound 1/2+ ground state 
and ^^Be a two-neutron separation energy higher than 1 
MeV. 

In scenario B, the energy of the 2si/2 state is given by 
the one-neutron separation energy in ^^Be. The energy 
of the lpi/2 state is fitted in order to reproduce at best 
the results for ^^Be and ^*Be. A good agreement with 
all quantities, including S'2n(^^Be) is found for a lpi/2 
state energy of 0.48 MeV, as shown in Table [Til This 
result for the lPi/2 state energy is in complete agreement 
with Refs.[l3, EOl- Moreover the 3si/2 state is found at 
1.33 MeV, corresponding to the 1/2+ state with £',.=:1.17 
MeV observed in Ref.fl^l. Note that we find for the 0+ 
state equivalent results as in Ref.[2| where the spin-orbit 
part of the DIS effective interaction was neglected. 

The pp-RPA amplitudes of 0+ state in ^'^Be obtained 
for scenarii A and B are presented in Table FVl Results de- 
pend strongly on the shell inversion hypothesis. As seen 
from Table [Vl in scenario A, Xajs amplitudes are small. 
This indicates that the core of ^^Be is little affected by 



TABLE V: Main pp-RPA amplitudes for Oj'" ground state in 
^*Be without (A) and with (B) inversion of 2si/2-lpi/2 shells. 



(2S1/2)' (ld5/2)^ 


Xq/3 
(lP3/2)^ (lPl/2)^ 


A -0.93 -0.49 


0.32 0.36 


(ld5/2)^ (lPl/2)^ (lPl/2 2pi/2) 


(IP3/2)' (2Si/2)2 


B -0.56 0.70 -0.63 


0.59 -0.45 



two-body correlations, contrary to results reported in Ta- 
ble |lll Now, looking at the amplitudes obtained for sce- 
nario B, a qualitative agreement with the components 
of ^^Be (see Table |n| is displayed. Indeed, Xap am- 
plitude for the configuration (2si/2)^ is quite large in- 
dicating that in ^^Be a configuration with two holes in 
the 2si/2-shell plays an important role. This result is in 
agreement with the amplitude for the (2si/2)^ configura- 
tion presented in Table |lll Thus concerning amplitudes 
for -'^'^Be ground state the scenario B gives solutions more 
consistent than scenario A. 




FIG. 3: Low lying spectra of Be obtained with pp-RPA 
without (A) and with (B) inversion in ^'^Be compared to ex- 
periment. 



In order to understand in more details consequences 
of scenarii A and B, we study excited states of ^°Be and 
"Be. 

In ^°Be, two 0+ and 1^ states are observed experimen- 
tally at E^(0^)=6.18 MeV [13 and E^(l^)=5.96 MeV 
[13, respectively. Two 2+ states at E^(2j^)=3.37 MeV 
and E2:(2j)=5.96 MeV [43] are also known from exper- 
iment. On the one hand, concerning scenario A, no 1^ 
state is predicted, whereas a Oj state at 6.31 MeV and 
a 2^ state at 3.83 MeV are present in the theoretical ex- 
cited spectrum. Even though those two states seem to 
be close to experiment, results may not have to be con- 
sidered satisfactory since they fail to obtain the lj~ state. 



On the other hand, scenario B displays a lj~ state at 3.82 
MeV and a Oj state at 3.96 MeV. The agreement is not 
quantitatively good. However the two states are close to 
each other as they are found experimentally. This sug- 
gests that these states are formed by two neutron holes 
coupled to the Oj excited state of ^^Be which has an ex- 
citation energy of 2.4 MeV; that is enough to shift the 
two states at the right energy. 

Results in ^''^Be from the two scenarii and compari- 
son with experiment are summarized in Fig|31 Here it is 
better suited to discuss results in terms of relative energy 
Er as in scenario A the experimental two-neutron separa- 
tion energy S2n is not well reproduced. In both scenarii, 
a good agreement for the relative energy E^ of the 0^ ex- 
cited state is obtained, Er=1.26 MeV and 1.24 MeV for 
scenario A and B respectively. The experimental value 
is equal to 1.22(18) MeV [l3|- The energies of the Oj 
state are close from each other in both scenarii. In Ta- 
ble |Vll the amplitudes for the 0^ state are displayed for 
both scenarii. In scenario A, O2 state is mainly built on 
(1^5/2)^, (2si/2, 3S1/2) and (2si/2)^ configurations. The 
values of Xap indicates the presence of correlations in 
the core nucleus. In scenario B, the Oj" state is explained 
with mainly (lpi/2)^ and (lpi/2, 2pi/2) configurations 
and the Xap are very small. Experimentally only the 
energy of this state is known. In a future experiment 
it would be interesting to investigate the spectroscopic 
factors of this state in order to discriminate between the 
two scenarii. 



TABLE VI: Main pp-RPA amplitudes for the QJ excited state 
in ^^Be without (A) and with (B) inversion of 2si/2-lpi/2 
shells. 



Xab 
(2Si/2)2 (2Si/2, 3Si/2) (ld5/2)' 


Xa;3 
(IP3/2)' (lPl/2)' 


A 0.42 0.58 -0.73 


0.27 0.27 


(lPl/2)^ (lPl/2, 2pi/2) 


(IP3/2)' (2Si/2)' 


B 0.73 0.64 


-0.08 0.04 



Concerning 1~ states in ^''Be, two low lying states with 
a relative energy £^j.=l-25 MeV and 1.32 MeV are found 
in scenario A. These values are a bit lower than the ex- 
perimental one Er=1.8±0.1 MeV obtained by Labiche et 
al. [4^. B{E1) transition probability between ground 
state and lj~ state has not been calculated as both states 
are not well reproduced in scenario A. 

Results appear to be much better in scenario B. In 
that case, the Ij^ state is found at Er=l.8 MeV {Ex^3.1 
MeV). Thus ground state properties as well as the energy 
of the lj~ state are very satisfactory. The El strength 
given by our calculation is B{E1)=3.7 xlO~'^ e^fm^. It is 
smaller than the strength obtained by Descouvemont et 
al. of 1.40±0.40 e^fm^ in their microscopic cluster model 
[43 . l48t . This difference comes from the fact that they 
do not assume any inversion in their model. They ob- 



tain a Si/2 state near threshold in their ^'^Be spectrum 
that enhances the El strength because of its spatial ex- 
tension. In FigHl the calculated El strength is shown 
for the different 1~ excited states obtained in scenario 
B. The main El strength is found in the energy-range 
Fjx^4:-7 MeV. The main strength is located at higher ex- 
citation energy than seen in ^^Bc. The same trend is 
observed by Sagawa et al. [23 ■ In Rcf.[40| Forsscn et 
al. have extracted the B(E1) distribution from the ex- 
perimental data of Ref ■ |48| . They have found a narrow 
B{E1) distribution peaked at about E^=2 MeV {Ex=5.3 
MeV). The shape of the strength is thus very different 
from the usual accumulation of El strength observed at 
low energy in other Borromean nuclei such as ^^Li and 
^He. Even using a phenomenological model of Coulomb 
dissociation with a lot of degrees of freedom, they did 
not manage to fit the shape of the El strength in ^"'Be. 
This result is a strong indication of a different structure 
in ^^Be in comparison with other Borromean nuclei. The 
absence of a low lying neutron s-state in the spectrum of 
^^Be and the appearance of a p-state is a possible expla- 
nation of the unusual shape of the soft El strength. In 
that sense results obtained for B{E1) arc in agreement 
with scenario B. 

Concerning the El transition probability strength, a 
value of 0.17 c^fm^ is found integrated below Er=3.2 
MeV {Ex=4:.5 MeV). This result seems low compared 
with the strength of 1.40±0.40 e^fm^ extracted by 



Forssen et al. Using the sum rule formula of Eq. (|16p 
with the calculated rms value of A, we obtain a value of 



Concerning the 2'^ state in ^^Be, it was first observed 
by Bohlen et al. [50| with an excitation energy of 
E:r=1.59(ll) MeV (E^=0.25(6) MeV). Then, this state 
was confirmed by Korsheninnikov et al. [35|. In a more 
recent experiment, Sugimoto et al. |ll| found the 2J 
state at E^,=l. 54(13) MeV (E,.=0.28±0.01 MeV). As al- 
ready discussed for ^^Be, this state is absent from our 
model as it is interpreted as an excitation of the ^^Be 
core. However, as shown in FiglS] a 2^ state at £^=1.6 
MeV is obtained in both scenarii. 

We think that further investigations are needed both 
on the theoretical and experimental parts in order to 
achieve a fully consistent description of ^^Be and ^^Be. 
In our study scenario B reproduces ground state as well 
as excited state properties, except the 2f state. These 
results are a strong indication on the validity of scenario 
B which assumes an inversion between 2si/2 and lpi/2 
shells in ^^Be. 



V. CONCLUSIONS 



We have employed the pp-RPA approach to describe 
even-even ®~^^Be isotopes from either a ^°Be or a ^^Be 
core. A WS potential corrected by a phenomenological 
particle-vibration coupling for the neutron-core interac- 
tion and the DIS Gogny force for the neutron-neutron 
interaction have been employed. 

Starting from the experimental spectrum of ^^Be and a 
^°Be core, our approach has been able to provide ground 



1.26 e^fm^ in agreement with the one of Ref.l4g|. The state properties as weU as excitation energies of the 0^ 



sum rule gives the El transition probability strength in- 
tegrated over the whole spectrum. We find that the El 
strength extracted below Er=3.2 MeV {Ex=A.5 MeV) is 
larger than the sum rule as found in Ref. [4^. Then the 
results deduced from the experiment of Labiche et al. 
may seem doubtful. 
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FIG. 4: Calculated El strength distribution in ^ Be for sce- 
nario B (with inversion). 



and lj~ states in ^^Be. As pp-RPA model assumes an 
inert core, it fails in describing the 2^ state that is prob- 
ably built on an excited ^°Be core. This issue could be 
cured in a further work introducing explicitly excitations 
of the core. The calculated B{E1) transition probability 
of the Ij^ state overestimates the experimental value. 

Concerning the description of the most exotic Beryl- 
lium isotopes, we have applied the same method with a 
^^Be core. In that case, ^"^Be spectrum has been needed 
as input of the calculation. Then, two scenarii have been 
tested: i) scenario A with a normal shell order in ^'^Be, ii) 
scenario B with a shell inversion similar to ^^Bc. Scenario 
A leads to several inconsistencies and is unable to repro- 
duce neither ground state nor excited state properties of 
-'^^Be. In scenario B, the low energy of the lpi/2 neutron 
state in the field of ^^Bc have been used as free parameter. 
To fix this parameter two constraints have been required. 
Firstly, the two-neutron separation energy in ^^Be should 
be in agreement with the experimental value. Secondly, 
^°Be and ^^Be ground and excited states should be in 
agreement with available data. These two different con- 
straints have led to the conclusion of a 2si/2-lpi/2 shell 
inversion in ^'^Be as observed in ^^Be and ^"Li. Such an 
assumption has implied the existence of a 1/2^ state in 
^■^Bc with an energy around 0.48 MeV, close to thresh- 
old. This state seems to have been observed in a recent 
experiment at RIKEN ^. 



The study of the two scenarii argues in favor of the 
conclusion of Refs.[l|, l^: the 2si/2-lpi/2 shell inversion 
is present in ^^Be. The present situation looks like the 
one of -'^^Li when theoretical studies on ^^Li had predicted 
the necessity to have an inversion in ^°Li [I3, l51[ , before 
it was confirmed experimentally [52, 153] . 



Appendix A: pp-RPA model and related algebra 

In order to simplify notations, we do not specify the 
total angular momentum in the pair creation and an- 
nihilation operators (see Eq.(I5|)), except when necessary. 
One obtains the following relations for amplitudes X and 
Y, 



Y.K?xr-Y.^Z'xi-p =w (Ai) 

ab ap 

ab a/3 

E^ir'<^-E^r^^i"^ =^-^- (A3) 

N M 



SaJfSX- (A4) 



N 



M 



We see from Eqs.® and (|5]) that the lowest eigenvalues 
among the N and AI solutions are directly related to 
the two-neutron separation energy in the {A + 2) and 
A nuclei, S2n{A + 2) and S2n{A) respectively. Indeed if 
Nq and Mo are the lowest energies of the two series of 
eigcnstatcs, then 



S2niA) = ftMa 
S2niA + 2) = -Qno- 



(A5) 

(A6) 



To determine the wave functions of the two systems, we 
introduce two operators Qjy and Q\j related to pair cre- 
ation and annihilation (see Eq.®). They are built in 
such a way that the core nucleus plays the role of a vac- 
uum. They verify the following equations. 



QjvlAO) = 
Qm\A,0) = 0. 



(A7) 
(A8) 



Moreover, wave functions of {A + 2) and {A — 2) systems 
are expressed as 



\A + 2,N) = qUA,6) 
\A~2,M) = QUA.O) 



(A9) 
(AlO) 



Requiring the orthonormalization of the wave function, 
and assuming the quasiboson approximation [17| , 

(A,0| [Aki,Al^] \A,0) « Sk,nSin (All) 

{A,0\[A,x,Al,]\A,b) « -S^^S^,, (A12) 



a relation between the operators Q (see Eqs. (|A9p and 
(|A10p ) and A (see Eq.(l5])) can be established. It reads 



Qi 



(^)/it 



k<l K,<\ 



K<A 

with the following properties 

Qn, QJv' 



k<l 

= Smm' 



(A13) 
(A14) 

(A15) 
(A16) 



All other commutators between two operators Q^ and Q 
are equal to zero. Inverting Eqs. (|A13[) and (|A14p . 



N AI 

^A-E^fQJv + E^ifs 



(A17) 
(A18) 



AT 



M 



Appendix B: Average value of one and two-body 
operators 



The calculation of the rms radius of Eq. (|T0|) needs the 
calculation of the average value of one- and two-body op- 
erators on the system formed by the two extra-nucleons. 
Then, given an operator F, the following evaluation is 
required 

(F) = {A + 2, Q\F\A + 2, 0) - {A, b\F\A, 0). (Bl) 

The wave function |v4 -f- 2, 0) corresponds to the ground 
state, 



1^ + 2,0) = Q5|A,0) 



(B2) 



SS-E^M.-E4°,Mi,, (B3) 

a<b a<P 



ith 



SO that Ea. (jBl[) can be transformed into 

(F) = {A,6\QoFQl\A,b) ' {ArO\F\ArO) 



= {A,b\Qo[F,Ql]\A,Q). 



(B4) 



The formula of Eq. (|B4[) can be applied to any one-body 
Fi or two-body F2 operators. For a one-body operator, 
one gets 



(F\) = 2^(m|Fi|n)^X(°)X 



(0) 
am 



2^(m|Fi|^)^xWxW, (B5) 



and for a two-body operator, one gets 



hi 



/C</,?71<TI 

2 Y. {kl\F2\jrv)Xf)xf^ 

k<l,fj.<L' 



(B6) 



K< A,/^<i^ 



Appendix C: Transition amplitudes in the A + 2 
nucleus 



In this Appendix, we consider a transition between the 
ground state (A^=0) and an excited state {N 7^ ) of the 
{A + 2) nucleus through a one-body operator (Fi), as for 
example electromagnetic ones. We do not specify angular 
momentum couplings for simplicity. The amplitude for 
such a transition is given by 

M{0 ^ N) =Y{^\Fi\j){A,0\QNala,Ql\A,0). (CI) 

The sum over i and j runs over all nucleon states, occu- 
pied or unoccupied. We can then rewrite A4 as 

M{0 ^ N) ^Y(t\F,\j){A,b\QN[ala,,Ql]\A,b). 



Then, using the following relation, 

(A6|SjvA1,-|A0> = {A + 2,N\Al\A,6) 



(C2) 



(C3) 



one obtains, 



X [{ai\Fi\ao)6bibo + {bi\Fi\bo)SaoaA 

Qo</3o,ai</3i 

X [(ai|Fi|ao)<5^o/3i + (/3i|Fi|/3o)i5aoai] 



(C4) 



In the particular case of an El transition between OJ 
ground state and the lj~ excited state and within the 
hypothesis of an inert core surrounded by two extra- 
neutrons, Fi is the electric dipole moment exciting the 
soft dipole mode [20| . 






(C5) 



i=l 



where e„ is the neutron effective charge defined in 
Eq.dini). Thus the soft El strength is 



-M2 



B{E1) = \M{Ot ^ ll)\ 



(C6) 
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